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ABSTRACT 

We report the discovery of a multiply lensed Lya emitter at z—3.90 behind the massive galaxy 
cluster WARPS J1415. 1+3612 at z=1.026. Images taken by the Hubble Space Telescope (HST) using 
ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet 
of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic 
observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide 
redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807 ± 
185 km s~ x . Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13 ± 
0.38") corresponds to a central velocity dispersion of 686T19 ^ m s_1 f° r the cluster, consistent with 
the value estimated from cluster member redshifts. Our mass profile estimate from combining strong 
lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results. 
Subject headings: 



1. INTRODUCTION 
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Surveys li ke Sp ARCS (jWilson et all &M i. SPT 
(|Ruhl et all 12004ft. the Red-Sequenc e Cluster Sur- 
veys C RCS. IGladderslj200l and RCS-2, ICladders fc Yeel 
[20051: lYee et al.l I2007T1 and the IRAC Shallow Survey 
( Eisenhardt et al.f2008l ) will discover thousands of galaxy 
clusters out to redshifts beyond z — 1. These large sur- 
veys will be used to probe dark energy and the underlying 
cosmology through direct measurements of the evolution 
of the cluster mass function. To fully utilize these data 
it will be essential to derive reliable mass estimates from 
the cluster properties measured in these surveys. 

Great effort has been dedicated to the development 
of cluster mass-observable relations in targeted cluster 
observations. Sunyaev-Zel'dovich (SZ) and X-ray ob- 
servations probe the hot ionized intracluster gas and 
have been used t o derive masses (|LaRoaue et alj|2006t 
lAllen et all 12008ft and scaling relations that tie mas s 
to the obse r ved S Z quantities (iBonamcn te et al.l 12008ft . 
iHicks et~a"l . (2006) derive scaling relations between ob- 
served X-ray properties and the optical richness of 
clusters while lEvrard et al.l (|2008ft model the relation- 
ship between velocity dispersion of cluster galaxies and 
cluster mass. Weak lensing measurements have been 
used to calibrate optical richness measurements from 
a large sample of clusters f r om the Sloan Digita l Sky 
Survey ([Koester et all l2007t Uohnston etHl 12007ft and 
to calibrate X-ray me asurement s through the mass- 
temperature relations jHoekstral 120071 : iBardeau et al.l 
l2007t iBerge et al.l 12008ft . Strong lensing measurements 
typically have much smaller errors compared with the 
other methods but are restricted to the central region 
of clusters and are rare. Nevertheless, strong lensing 
systems have been studied in small samples of clusters 
(jComerford et al . 1 12001 [Hennawi et al.ll2008ft . 

Large cluster surveys out to z ~ 1 have the poten- 
tial to constrain the dark energy equation of state (e.g. 



2 



lVoitJ |2"005). provided that the mass-observable relations 
are well calibrated. Full exploitation of the more distant 
clusters requires deep space-based observations. As a re- 
sult, mass-observable relations have been studied in far 
fewer clusters at z > 1 relative to low redshift clusters. In 
a 219 orbit program (Program number 10496, PI: Perl- 
mutter) to search for supernovae (SNe) with the Hubble 
Space Telescope (HST), we observed 25 of the highest 
redsh ift galaxy clusters known at the time (jDawson et al.1 
2009). These images support a rich program of clus- 
ter studies including the calibration of mass proxies at 
z > 1. One particular cluster from this program, se- 
lected from the Wide Angle ROSAT P ointed Survey, 
WAR PS J1415. 1+3612 at z = 1.026 (iPerlman et al. 
2002), has already been studied in X-ray dj Vlaugha n et al.l 
12009: lAllen et all l200l ) and SZ (jMuchovej et all l2007h 
observations. We now have deep images from the Ad- 
vance Camera for Surveys (ACS) and spectroscopy from 
the Faint Objec t Cam era and Spectrograph (FOCAS: 
iKashikawa et al.l [2002) on Subaru that eveal a pro- 
nounced strong lensing arc of a source Lya emitting 
galaxy at z = 3.90. The new data from this program 
combined with previous and ongoing measurements en- 
able a multi-probe analysis of the cluster mass-observable 
relation for this high redshift cluster. 

Here we present an analysis of the strong lensing and 
dynamical mass of WARPS J1415. 1+3612. The letter is 
organized as follows: we describe the observations and 
data in £j2j mass estimates are derived and compared in 
iJHl and the summary is found in f]4] Throughout this 
letter we use AB magnitudes and assume a cosmology 
with fl M = 0.3, fi A = 0.7, and h = 0.7. 

2. OBSERVATIONS 
2.1. ACS Images 

WARPS J1415. 1+3612 was observed 7 times with ACS 
from November 2005 through April 2006 for a total inte- 
gration of 2425 seconds in the F775W filter and 9920 
seconds in the F850LP filter, hereafter 1775 and zsso- 
Individual exposures w ere coadded using MultiDrizzle 
(|Fruchter fc Hookl l2002T i at a resolution of 0.05"/pixel. 
We use 25.67 8 and 24.867 as zero points for 1775 and z%5o 
respectively ()Sirianni et al .1120051 ). 

The deep ACS images revealed a complex strong lens- 
ing system near the cluster core. The composite color 
image from the 1775 and zsso data is shown in Figure [TJ 
The arc system consists of at least five images. The 
main arc (A) is 6.75" (54.4 kpc at the cluster redshift) 
from the center of the cluster, which we take to be the 
position of the brightest cluster galaxy (BCG). To the 
SW of the main arc, there is a triplet of arcs (B, C and 
D) around a spectroscopically confirmed cluster member 
(#13). About 5" south of the triplet, there is a fifth 
arc (E). The total 1775 isophotal magnitudes for the arcs 
A - E are 23.44, 25.75, 25.46, 25.13 and 25.62, respec- 
tively. We use the 4775 magnitudes because they suffer 
less contamination than the z^^q magnitudes from the 
much redder elliptical cluster members. The magnitudes 
for the arcs B, C and D are obtained after the lensing 
galaxy has b e en su btracted using the b-spline model of 
iBolton et all ((2001 . These five images lie very close to 
an imaginary arc centered on the BCG that subtends 
slightly more than 90° . The radius of the Einstein ring 



for the cluster potential is taken to be the average of the 
two circles mentioned in Figure[IJ 7.13 ± 0.38". The un- 
certainty in the radius is determined from the difference 
between the two circles and the average. 

The B band luminosity for this cluster is estimated by 
summing up light from all the galaxies within approxi- 
mately 1 Mpc of the BCG and then su btracting a back- 
groun d we estimate from the GOODS (|Giavalisco et al.1 
l200l images. We apply a i-T-correction to transform 
Z850 to restframe B magnitude. The total luminosity is 
determined to be L B = 2.92 ± 0.88 x 10 12 L Bo . 
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Fig. 1.— Top: Central Region of WARPS J1415. 1+3612. Strong 
lensing arcs are labeled "A"-"E" and cluster galaxies of interest are 
labeled "#13" and "#15". The radius of the Einstein ring is taken 
to be the average of the radii of two imaginary circles (not shown), 
both centered on the BCG. The first circle with radius 6.75" passes 
through the center of arc A and the second, with radius 7.51", 
passes through cluster member #13 (the lensing galaxy). F and G 
are discussed in the text. Bottom: The panel on the left shows 
smoothed spectra of arc A; the triplet B, C and D; and arc E. Zero 
flux is shown as the horizontal dotted lines. The diamonds with 
error bars are the average values of flux in 500 A bins. The differ- 
ences between the averages of the flux for wavelengths longer and 
shorter than the peak for arcs A and E (dashed lines) indicates the 
detection of discontinuity around the emission feature. The panel 
on the right shows the same spectra smoothed to the spectrograph 
resolution (FWHM=600 km s -1 ) and resampled by 200 km s — 1 
pixel in the log linear scale. 
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2.2. Subaru Spectroscopy 

Spectroscopic observations were made with the FO- 
CAS spectrograph on the Subaru telescope. SN candi- 
dates provided the primary targets for spectroscopic ob- 
servations. The strong lensing arcs described in H2.1l and 
galaxies with 1775-2:850 color consistent with early type 
galaxies at the cluster redshift provided secondary tar- 
gets for multi-slit spectroscopy. We obtained spectra for 
a total of 138 objects and measured redshifts for 95 of 
them. Six masks were used, four with the 300R/SO58 
grism and filter and two with the 300B/SY47 grism and 
filter. The total integration time is approximately two 
hours for each mask. Arcs B, C and D were observed 
with the former setup, resulting in spectra that start at 
5800A, while arc E was observed with the latter, result- 
ing in spectra that cover a bluer wavelength range. Arc 
A was observed with both setups. The slit width is 0.8", 
corresponding to a spectral resolution of R ~ 500. Fur- 
ther details can be found in Morokuma et al. 2009 (in 
preparation) . 

Redshifts for cluster members are obtained by fitting 
a series of galaxy spectral eigenfunctions to the spectra 
weighted by the error spectra. For the majority of the 
spectra, the 3933 A and 3968 A Ca II H and K absorption 
lines are the dominant feature detected, along with the 
4000 A break and the [Oil] emission line. Each redshift 
is assigned a quality flag Q, with Q = 3 for those with 
multiple strong features, Q = 2 for those with a single 
strong feature and Q = 1 for those with only weakly de- 
tected features. Cluster members selected from the spec- 
troscopically confirmed galaxies inside a radius of ~ 2.1' 
are listed in Table 1. The redshift of the cluster is deter- 
mined to be 1.026 from the average of the cluster member 
redshifts listed in the table. At this redshift, a radius of 
2.1' corresponds to approximatly 1 Mpc. We calculate 
the cluster velocity dispersion from the 21 g alaxies with 
Q > 1 . Using the biweight scale estimator of lBeers et al.l 
(1990), we estimate er v = 807 ± 185 km s" 1 in the cluster 
rest frame. The measurement uncertainty is determined 
from bootstrap resampling of the data. 

Spectra from each of the arc candidates (A - E) show 
significant emission around 5950 A (Figure Q}. The con- 
sistency in the observed emission wavelengths, surface 
brightness, color, and morphology of these five objects, 
and their arrangement around a circle centered on the 
BCG all strongly suggest that this is a multiply imaged 
system at z > 1.026. In this case, the line is almost 
certainly Lya at restframe 1215.67 A. The detection of 
discontinuity in flux for wavelengths bluer and redder 
than the emission feature in the spectra for arcs A and E 
lends support to the identification of this line with Lya. 
Assuming this emission is due to Lya, the redshifts of 
each arc are reported in Table [1] As indicated by the 
apparent discrepancy in redshifts, there is ~ 9 A (or 
453 km s _1 in the rest frame of the cluster, Figure [l| 
difference between the observed wavelengths of Lya in 
arc A and in the triplet (arcs B, C, and D). This is 3 
times the statistical error, which takes into account the 
spectral resolution. This difference is likely due to the 
placement of the arcs B, C and D near the edges of the 
slit and contamination from the lensing galaxy (cluster 
member #13). 



TABLE 1 
Spectroscopic redshifts 



Object ID 


R.A. (J2000) 


Dec. (J2000) 


z 


Q 


1 


14:15:13.08 


+36:12:25.4 


1.0262 ± 0.0003 


3 


2 


14:15:12.68 


+36:12:59.6 


1.0331 ± 0.0003 


3 


3 


14:15:11.27 


+36:14:35.7 


1.033 ± 0.001 


2 


1 


14:15:11.64 


+36:11:37.8 


1.028 ± 0.002 


1 


5 


14:15:07.08 


+36:14:15.6 


1.0255 ± 0.0003 


3 


6 


14:15:07.68 


+36:13:06.3 


1.025 ± 0.001 


2 


7 


14:15:13.95 


+36:12:37.8 


1.017 ± 0.001 


2 


8 


14:15:10.36 


+36:11:31.2 


1.0256 ± 0.0003 


3 


9 


14:15:15.91 


+36:12:58.9 


1.0187 ± 0.0003 


3 


10 


14:15:10.66 


+36:11:39.8 


1.028 ± 0.001 


2 


11 


14:15:08.08 


+36:12:21.0 


1.0234 ± 0.0003 


3 


12 


14:15:10.54 


+36:12:09.9 


1.0264 ± 0.0003 


3 


13° 


14:15:10.59 


+36:12:07.9 


1.023 ± 0.002 


1 


1 1 


14:15:10.19 


+36:11:49.6 


1.0313 ± 0.0003 


3 


15" 


14:15:11.12 


+36:12:04.0 


1.0252 ± 0.0003 


3 


16 


14:15:11.47 


+36:12:22.8 


1.018 ± 0.001 


2 


17 


14:15:06.61 


+36:12:18.0 


1.033 ± 0.001 


2 


18 


14:15:09.12 


+36:13:09.1 


1.0260 ± 0.0003 


3 


19 


14:15:13.45 


+36:12:56.9 


1.035 ± 0.002 


1 


20 


14:15:12.60 


+36:11:09.7 


1.023 ± 0.001 


2 


21 


14:15:00.02 


+36:13:02.8 


1.011 ± 0.001 


2 


22 


14:15:05.89 


+36:13:03.1 


1.025 ± 0.001 


2 


23 


14:15:13.00 


+36:13:07.0 


1.0301 ± 0.0003 


3 


24 


14:15:08.41 


+36:13:04.6 


1.035 ± 0.002 


1 


25 


14:15:10.82 


+36:10:33.0 


1.023 ± 0.001 


2 


A 


14:15:10.80 


+36:12:09.1 


3.8983 ± 0.0005 


2 


B-D 6 


14:15:10.57 


+36:12:06.9 


3.891 ± 0.002 


2 


E 


14:15:10.58 


+36:12:00.8 


3.897 ± 0.001 


2 



a Highlighted in Figure 1 

Triplet unresolved in 2-d spectra. 



Preliminary spectroscopic results show that object F 
in Figure [T] is a possible counter arc of the same source 
galaxy. Object G in FigureQ]has the same color as arcs A 
- E and lies near the Einstein ring but spectrocopic data 
are inconclusive as to whether it has the same redshift 
as the confirmed arcs. 

3. DISCUSSION 

The total mass inside the Einstein ring of 
r E = 57.5 ± 3.1 kpc is 1.96±g;jg x 10 13 M o (e.g. 
INaravan fe BarteTm ann 1993). If we assume a singular 
isothermal sphere (SIS) profile for the central region 
of the cluster, this mass corresponds to a rest frame 
velocity dispersion of 686l 19 km s^ 1 , consistent with 
the value estimated from cluster member redshifts. We 
can estimate the magnification f actor for arc A by using 
H = 1 1/(1 — rg/r) | (e.g. INaravan fc Bartelmannl [19961) 
for the SIS profile, where r is the distance of the image 
from the lens center. If we take te = 7.13" (or 57.5 kpc), 
/i ~ 18 for arc A. This implies that the unmagnified 1775 
magnitude for arc A is ~ 26.6. 

To get a qualitative estimate of the change in the de- 
rived strong lensing mass caused by relaxing our assump- 
tion of a spherically symmetri c mass distribution, w e con- 
struct elliptical NFW models (|Glose fc Kneibll2002ft with 
a critical curve that passes through at least part of arcs 
A and E and passes between arc C and the lensing galaxy 
(cluster member #13). The model with the largest ellip- 
ticity that still meets these requirements has an elliptic- 
ity of 0.1. The mass enclosed within its critical curve is 
1.79 x 10 13 M Q , within the errors of the spherical models. 

The dynamical mass is estimated from the redshifts 
of the 21 cluster members with Q > 1 in Table [T] us- 
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ing the viria l scali ng relation from the simulations of 
lEvrard et alj 



&dm(M,z) = CTDM,1E 



h(z)M 2 , 
10 15 Af ( 







(1) 



where auM(M, z) is the dark matter velocity dispersion, 
h(z) = H(z)/100 km s _1 , and M 2 oo is the mass within 
^200 (the radius of a sphere where the mean density is 
200 times the critical density at the redshift of the clus- 
ter). The constants <7dm,15 = 1082.9 ± 4.0 km s _1 and 
a = 0.3361 ± 0.0026 ar e determined from the simulations 
of lEvrard et al.l (|2008f ). Assuming that the velocity bias 



b v = crv,gai/crDM = 1, we find M 200 = 3.33+tiB x 1O 14 M 
and r2oo = 0.97 ±0.22 Mpc. Thus the mass to light ratio 
is ~ 114 Mp,/L BP) . In Figure El we fit an NFW profile 
([Navarro. Frenk. fc Whital997F to the total cluster mass 
using the strong lensing and dynamical mass estimates 
The model has a concentration parameter c 
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and a core radius r s = 0.179ig'o55 Mpc, where uncer- 
tainties are computed via Monte Carlo simulation. 

For a simple test of the mass-observable relation in 
WARPS J1415.1+36, we compare our strong lensing and 
dynamical estimates and a model of the mass profile to 
estimates of weak lensing mass, X-ray mass, and SZ mass 
in Figured] Using the same deep ACS images, Jee et al. 
(2009, in preparation) derive the weak lensing mass as- 
suming a n NFW profile with a mass-concentration rela- 
tion from |Bu11ocOOI3 (HQQB)- The X-ray mass M 2 5 00 ( z ) 
(|Maughan et al.ll2006f ) is determined from a best fit two- 
dimensional elliptical /3-model analyzed at r 2 5oo(z) ■ The 
SZ value of M 2 5oo(z) is derived from a spherical, isother- 
mal /3-model fit to 3 GHz data from the SZ Array 
(|Muchovei et alj [20071 ). In the X-ray and the SZ mass 
estimate, A(z) is a r edshift dependent dens ity contrast 
as described in e.g. iMaughan et alj ([2006D . Both the 
weak lensing mass and the X-ray mass are consistent 
with our results. However the SZ mass is lower than the 
prediction of our mod el by roughly a factor of two. In 
iMuchovei et all (|2007f) . it is noted that the SZ mass es- 
timate is marginally lower than the X-ray value, likely 
due to low S/N in the SZ measurements, uncertainties in 
the absolute calibration, and possible systematic errors 
from the model assumptions. Additional observations 
with the SZA at 30 GHz and 90 GHz should improve the 
constraints on the SZ mass estimate and shed light on 
the apparent discrepancy. 

Using numerical simulations, iDolag et al.1 (|2004f ) pre- 
dict the dependence of the average halo concentration 
parameter on the cluster redshift and virial mass 



c(M,z) 



Co 



1 + z 



M 



(2) 



I Shaw et al.l (|2006f ) find the best fit parameters to be 
c = 6.47 ± 0.03 and 7 = -0.12 ± 0.03 from N-body 
simulations. Here M is the v irial mass and M* = 
11.0 x 10 14 M Q ([Shaw et al.ll200l . If we assume the dy- 
namical mass (Af2oo=3.33 x 10 14 M Q ) to be the virial 
mass, Equation [2] predicts the average concentration to 
be c = 3.70lQ 34ti.3, where the first set of errors are 
from the errors in the best fit parameters (co and 7) 
and the second set is from the scatter about the average 



( Bull ock et al.l [200l[ ) . Our concentration of c = 4.96 is 
higher than the predicted value of c but is well within its 
errors. 

We now examine the strong lensing arcs B, C, and 
D around cluster member #13. The radius of the Ein- 
stein ring for the lensing galaxy is estimated by the dis- 



tance between the galaxy and arc C, 0.65 



//+0.35" 

-0.30" 



The 



upper and lower bounds are set by the distances from 
the lensing galaxy to arcs D and B, respectively. Assum- 
ing that the deflection angle due to the cluster potential 
does not vary appreciably for arcs B, C and D, the mass 
enclosed within the Einstein radius is estimated to be 
1.65l?;?7 x 1011 M o- 
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Fig. 2. — Comparison of different mass estimates for WARPS 
J1415+36. The strong lensing, weak lensing and dynamical mass 
estimates presented in this work are labeled above. M 2 5oo( z ) from 
the X -ray UMauehan et al.1 |200D and SZ data [IMuchovei et al.1 
2007) are also shown. The dashed line is an NFW profile fit to the 
strong lensing mass at rg=0.0575 Mpc and the dynamical mass at 



r 20 o=0.97 Mpc (c = 4.96; 



- 1 03 



: 0.179 



+0.078 



Mpc). 



4. SUMMARY 

The galaxy cluster WARPS J1415+36 at zi ens = 
1.026 is o ne of the most distant s trong lensing clusters 
reported ([Thompson et~aH I200U IGladders et all 120031 : 
Gilba nk et alj 120081) . Although the SZ mass estimate 
is not consistent with the rest of the data for this clus- 
ter, this work suggests that the mass-observable relations 
derived from numerical simulations and observations of 
strong lensing, weak lensing, dynamical mass, and X-ray 
mass still h old for clusters a t reds hifts at z ~ 1. 

Recently Gilba nk et alj (|2008f ) report multi-probe 
mass estimates for the strong lensing cluster, RCS2319, 
at z = 0.91. Future surveys will soon reveal many new 
high redshift, strong lensing systems and provide deep, 
multiwavelength images of these clusters. With the cor- 
responding improvements in mass-observable relations in 
the decelerating regime, galaxy clusters will soon become 
an even more effective tool for cosmology. 
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